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ABSTRACT 

Aims. We identify and characterize low-mass stars in the ~ 3 Myr old Trumpler 16 region by means of a deep Chandra X-ray observation, and 

study their optical and near-IR properties. We compare X-ray activity of Trumpler 16 stars with known characteristics of Orion and Cygnus OB2 

stars. 

Methods. We analyzed a 88.4 ksec Chandra ACIS-I observation pointed at the center of Trumpler 16. Because of diffuse X-ray emission, 

source detection was performed using the PWDetect code for two different energy ranges: 0.5-8.0 keV and 0.9-8.0 keV. Results were merged 

into a single final list. We positionally correlate X-ray sources with optical and 2MASS catalogues. Source events were extracted with the IDL- 

based routine ACIS-Extract. X-ray variability was characterized using the Kolmogorov-Smirnov test and spectra were fitted by using XSPEC. 

X-ray spectra of early-type, massive stars were analyzed individually. 

Results. Our list of X-ray sources consists of 1035 entries, 660 of which have near-IR counterparts and are probably associated with 

Trumpler 16 members. From near-IR color-color and color-magnitudes diagrams we compute individual masses of stars and their A, values. 

The cluster median extinction is Ay =3.6 mag, while OB-type stars appear less absorbed, having Ay =2.0 mag. About 15% of the near-IR 

counterparts show disk-induced excesses. X-ray variability is found in 77 sources, and typical X-ray spectral parameters are Nh ~5. 37x10"' 

cm"^ and kT~1.95 keV, with Icr dispersions of 0.45 dex and 0.8 keV, respectively. OB stars appear, softer with a median kT~0.65 keV. The 

median X-ray luminosity is 6.3x 10^^" erg s"' , while variable sources show a larger median L^ value of 13x 10^^" erg s"' . OB-stars have an even 

higher median Lx of 80x10^° ergs"', about 10 times that of the low-mass stars. 

Conclusions. The Trumpler 16 region has a very rich population of low-mass X-ray emitting stars. An important fraction of its circumstellar 

disks survive the intense radiation field of its massive stars. Stars with masses 1.5-2.5 Mq display X-ray activity similar to that of stars in 

Cyg OB2 but much less intense than observed for Orion Nebula Cluster members. 
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1. Introduction 

The Carina nebula region (NGC 3372) is one of the most mas- 
sive star formation regions of the Galaxy. It is associated with a 
giant Hii region spanning about 4 deg^ of the sky and being bi- 
sected by a prominent V-shaped dark gas and dusty lane. This 
prominent young structure is not as compact as some of the 
other young galactic clusters, but seemingly related to a spiral 
feature. In this direction, we are looking almost tangentially to 
the now recognized Carina-Sagittarius spiral arm, at the edge of 
a giant molecular cloud extending over about 130 pc and with 
a content in excess of 5x10^ solar masses dGrabelskv et al 



19881) . The concentration of massive stars (i.e. M>20 Mq) in- 
teracts with the parent giant molecular cloud of the region, 
leading to trigger ed star formation e vents on intermediate to 
lower masses (e.g. lSmith et al.ll2004i) . 
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This region harbors several open clusters and/or star concen- 
trations (Trumpler 14, 15 and 16; Collinder 228 and 232; 
Bochum 10 and 11) c ontaining more than 60 known O-type 
stars (Feinsteiii, Il995h . Large cavities within the giant molec- 
ular cloud are supposed to be carved out by the Tr 14 and 
16 open clusters, which contain most of massive stars of the 
region. In particular Tr 16 includes three rare main-sequence 
03 stars, the Wolf-Rayet (WR) star HD93162 and the famous 
luminous blue variable (LBV) //Carinae . There is a histori- 
cal controversy about the distance and age of Tr 14 and Tr 16 
(IWalbornL Il995[). For instance, from ext ensive spectroscopy 
and photometrv lMassev & JohnsonI (Il993h find 3.2 kpc for both 
clusters. However, photometric studies are strongly affected 
by differential extinction in the region and peculiar redden- 
ing, and_so_tfie derived distance are different. An example is 
the ICarraro et a l. (2004) work, who for different R=Av /E(B- 
V) values (3.48 and 4. 16 for Tr 16 and Tr 14 regions), compute 
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distances of 4.0 kpc and 2.5 kpc, respectively. A more reliable 
distance (2250+180 pc), was derived from proper motion and 
Doppler velocities of the expanding ;; Carinae Homun culus us- 
ing HST-STi:^ observations (David son & Humphrey s. 1997). 
Recent work dTapia et all 120031) derives a common distance 
DM=12.14 (2.7 kpc) and an age between ~1 Myr and 3Myr, 
for Trl4 and Trl6, respectively. For this study, we adopt for 
Trumpler 16 a distance of 22 50 pc and a n age o f 3 Myr. This 
young age agrees with the ISmith et al.l (120001) results, who 
report the existence of several emb edded IR sourc es where 
star fo rmation might be active. Also. iDeGioia-Eastwood et al.l 
(120011) confirm clear evi dence of pre-main s equence (PMS) 
stars in the region, while Brooks et al.l (1200 ih have identified 
two compact Hn regions possibly linked to very young O-type 
stars. Finally. iHagele et al.l (120041) report a compact cluster of 
infrared PMS-stars in Tr 16. 

Of the existing methods to identify young stellar popu- 
lation s, the use of X -ray emission is perhaps the least bi- 
ased (IFeigelson et al.l 2002). While in main-sequence (MS) 
stars, from late A to M dwarfs. X-rays are believed to orig- 
inate from the hot coronal gas that is heated by stellar dy- 
namo magnetic fields (Maggio et al., 1987), for late type Pre- 
MS stars (T-Tauri stars (TTSs)) X-ray emission is attributed 
to solar-like coronal activity but elevated by a factor of 10"^- 
10^ (Feigelson & Montmerle, 1999). Several authors suggested 
the possibility of detecting early pre-main sequence (PMS) ob- 
jects through their hard X-ray emission escaping the highly 
obscured regions (see Walt e r (1992); Ka r nata et al. (1997); 
Hofner & Churchwelll (Il997l) ; iHofner et all (|2002|) ). Recendy, 
X-ray surveys have been successful in identifying the young 
and pre-MS population in star-forming regions, including: i- 
deeply embedded Class 1 young stellar objects (YSOs), ii- 
low-mass T-Tauri PMS stars. Hi- intermediate-mass Herbig 
Ae/Be PMS stars, iv- zero-age MS stars. Moreover, X-ray 
emission from low-mass pre-MS stars usually exhibits a strong 
variability that helps to confirm membership. 

On last decade. X-ray observations of young stars on star- 
forming regions were intensified thanks to the high spatial 
resolution and the improved broad-band ([0.2-12.0] and [0.5- 
10.0] keV) effective area of the XMM-Newton and Chandra 
satellites. A first X-ra y surv ey in the Carina region by 
(l2003h was performed on the basis 
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of two early XMM-Newton observations (rev#l 15 and #116) 
centered on ;; Carinae . Because of the spatial resolution of the 
EPlCo camera and relatively short exposure time of the obser- 
vations (~35 ksec), they detected only 80 X-ray sources, most 
of them related to the massive OB-type stars with L^ ~10^^- 
lO^'* erg s '. Before the observation used here, three Chandra 
observations were obtained on this region, two (obsld50 and 
1249) in the timed exposure mode, and the third (obsId51) 
in the continuous clocking mode, which produces no image. 
Using only observation obsld 1249, Evans et al. (2003) pre- 
sented luminosities and hardness ratios of the hot stars in Tr 16, 



' Data from Hubble Space Telescope (HST) with the Space 
Telescope Imaging Spectrograph (STIS) 

^ European Photon Image Camera has about six times less spatial 
resolution than Chandra ACIS-I camera. 



and part of Tr 14. Low -re s olution X-ray spect ra of luminous 
sources were discussed by Evans et al.l ( 120041) . However the 
short exposure time of such observation (~ 9.5 ksec) was a se- 
rious limitation for the study of intermediate- and low mass 
stellar population of the region. This limita tion exists even if 
obsld. 50 and obsld. 1249 are combined (iSanchawala et al.L 
2007), reaching completeness just at X-ray luminosity (Lx ) of 
~7xl0"^' erg s"', i.e. the X-ray emission level typical of single 
O- and early B-type stars. 

In this paper we present results of the analysis of the deep- 
est X-ray observation ever done in this region (~90 ksec). 
Section |2] gives details on the observation and data reduction 
procedures. Section 3 explains the method used to detect the 
sources, photon extraction and the construction of the catalog. 
In section 4 we present results of the cross-correlation with ex- 
isting near-lR and optical catalogs of objects and their charac- 
terization based on their color-color (CC) and color-magnitude 
(CM) diagrams. Section 5 presents a statistical study of vari- 
ability in the X-ray domain. Section 6 is dealing with results 
of the analysis of extracted X-ray spectra. In section 7 we dis- 
cuss X-ray luminosities of stars and compare them statistically 
with the X-ray source population of ONC and Cygnus 0B2 star 
forming regions. In section 8 we discuss X-ray and stellar pa- 
rameter of O- and early B- type stars. Finally, in section 9 we 
give a summary and draw conclusions of the paper. 



2. X-ray observations 

Trumpler 16 was observed with the ACIS detector on board the 

Ion 



Chandra X-ray Observatory (CXO) (IWeisskopf et al.[ l2002h 



2006 August 3 l|j (obsld 6402), as part of the Guaranteed Time 
Observation (GTO) Chandra program. The total effective ex- 
posure time was 88.4 ksec. The data were acquired in very 
FAINT mode, to ease filtering of non-X-ray events, with six CCD 
turned on, the four comprising the AClS-1 array [0,1,2,3], plus 
CCDs 6 and 7, part of ACIS-S. However, data from the lat- 
ter two CCDs will not be used in the following because of the 
much degraded point spread function (PSF) and reduced effec- 
tive area resulting from their large distance from the optical 
axis. The AClS-1 17'xl7' field of view (FOV) is covered by 4 
chips each with 1024x1024 pixels (scale 0.49" px"'). The ob- 
servation was pointed toward R. A.= lO'' 44™ 47.93' and DEC=- 
59° 43' 54.21", chosen to maximize the number of stars in the 
FOV and close to the optical axis, but also including most of 
the OB stars of the cluster Figure [T]-left shows Trumpler 16 as 
seen in X-rays by our ACIS-I observation. 

Most of the observed X-ray sources in Fig[T]-left are likely 
located toward the central part of Trumpler 16 , spatially con- 
strained by the dark V-shaped dust lane of the Carina region 



(IBrooks et al.L 119981) . However, note the small X-ray source 
concentration inside this dark structure, located towards the 
South-East part of Fig 1 . Deep near-IR observation are needed 
to reveal counterparts of such population. We defer to a forth- 
coming paper the use of some public HST - WEPCI^J obser- 



^^ Observation start date is JD 273350372.4528. 

* Hubble Space Telescope - Wide Filed Planetary Camera 2 
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Fig. 1. Left: Color-coded ACIS-I image of the 17' x 17' field in Trumpler 16 (see color version in the electronic edition). Kernel 
smoothing was applied to highlight point sources. Energy bands for the RGB image are [0.5:1.5], [1.5:2.2], and [2.2:8.0] keV for 
the red, green, and blue colors, respectively. Contour lines show the spatial distribution of the dark V-shaped dust lane. Right: 
ACIS-I image composed with a Hq, image. This image clearly show the sharp cut, along the V-shape structure, in the spatial 
distribution of X-ray sources. 



vations to find signatures of star-environment interactions for 
some of the X-ray sources. 

2.1. Data reduction 

Data reduction, starting with the Level 1 event list provided by 
the pipeline processing at the CXO, was performed using CIAO 
3.3.0.10 and the CALDB 3.1.0 set of calibration files. We pro- 
duced a level 2 event file using the acis_process_event CIAO 
task, taking advantage of the VF-mode enhanced background 
filtering, and retaining only events with grades = 0,2,3,4,6 and 
status=0. Photon energies were coiTected for the time depen- 
dence of the energy gain using the corr_tgain CIAO task. 
Intervals of background flaring were searched for, but none 
were found. We will hereafter assume a non-variable back- 
ground. To improve the sensitivity to faint sources, given the 
spectrum of the background and that of typical sources, we fil- 
tered out events outside the [500:8000] eV energy band. 

3. Data analysis 

We built images in the three band-passes 0.5-1.5 keV (soft), 
1.5-2.2 keV (medium) and 2.2-8.0 keV (hard). Before the 
color-coded image combination, we corrected single band- 
images by variations in exposure, sensitivity and vignetting, by 
computing and dividing by its respective exposure maps. We 
construct a color-coded X-ray image of the region by compo- 
sition of the three soft (red), medium (green) and hard (blue) 
images (see Fig.[T]l- This image is a 17xl7arcmin field which 



comprises the center of Tr 16 and marginally the southeast part 
of Tr 14. In addition to revealing a huge number of X-ray point- 
like sources with different color-energies, inspection of Fig.[I] 
suggests the preser ice of soft diffuse X-ray emission in the re- 
gion. In particular. iTownsIevi (120061) use a 57 ksec ACIS-I ob- 
servation centered on Tr 14 to study the diffuse X-ray emis- 
sion in such a region and in the North- West part of Tr 16. They 
explain the diffuse X-ray emission as likely arising from the 
fast O-star winds that shock and thermalize the surrounding 
medium. Because homogeneous data analysis techniques are 
required to do justice to a comparison of diffuse X-ray emission 
in this data-set with that of Tr 14, this study will be presented 
in a forthcoming paper. 

3.1. X-ray source detection 

Source detection was perfor med with the Palermo Wavelet 
Detection code, PWDetec|l (JDamiani et aI.L Il997bl) . It ana- 
lyzes the data at different spatial scales, allowing the detection 
of both point-like and moderately extended sources, and effi- 
ciently resolving close source pairs. The most important input 
parameter is the detection threshold (SNR), which we establish 
from the relationship between background level of the observa- 
tion and expected number of spurious detections due to Poisson 

noisqj- 

Diffuse soft (~0.2-1.2 k eV) X-ray emission has already 
been identified in this region JCorcoran et aUll995tlTownslev , 



http://cxc.harvard.edu/ciao/ 



^ See http://www.astropa.unipa.it/progetti_ricerca/PWDetect 
' This last quantity was determined from extensive simulations of 
source-free fields (see. lPamiani et aUll997al) 
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20061) and causes different (non-uniform) X-ray background 
levels across the FOV of our observation. The best way to re- 
duce the diffuse emission contribution is to discard soft photons 
in the detection procedure. We compute background (BKG) 
levels in four different energy ranges: 0.5-8.0 keV, 0.7-8.0 keV, 
0.9-8.0 keV and 1.1-8.0 keV. Since exposure maps are needed 
by the source detection algorithm, we used event files in these 
energy bands as input to the ClAO tool mkexpmap and assumed 
a monochromatic spectrum (kT=2.0 keVn The background 
level for each energy band was computed as the average of val- 
ues measured in three circular regions free of sources. If we ac- 
cept 10 spurious detections in the ACIS FOVq mean computed 
background levels can be translated into different SNR thresh- 
olds for each energy band, i.e.: 4.7, 4.65, 4.6 and 4.55, respec- 
tively. These input parameters for PWDetect result in different 
numbers of detected sources in each energy band, for instance, 
1266 sources in the 0.9-8.0 keV band, greater than detected in 
the "canonical" 0.5-8.0 keV band (1214 sources). This is con- 
sistent with the effect of diffuse soft X-ray emission in masking 
weak sources. 

A careful visual inspection was performed on the 0.5-8.0 
keV and 0.9-8.0 keV source lists. We manually rejected re- 
spectively 242 and 271 detections, considered spurious either 
because they were produced by different instrumental artifacts 
(e.g. CCD gaps, detector edges, false detections along the read 
out trails), or since they resulted from the "fragmentation" into 
discrete sources of the extended emission making a toroidal 
ring around the LBV star ?7Carinae (Fig. [TJ. Furthermore ex- 
amination of afterglowl"l contamination led us to discard 29 
and 20 false detections, respectively. The final source list was 
constructed by merging both catalogues: the list of 943 sources 
detected in the 0.5-8.0 keV band and the 1004 sources de- 
tected in the 0.9-8.0 keV band. Merging was performed us- 
ing a criterion of maximum detection significance. The two 
catalogues coincide in 798 sources, consisting of 322 sources 
with Sigo.5-8.o>Sigo.9-8.o, while for the remaining 476 sources 
Sigo.9-8.o>Sigo.5-8.o (where Sigj,--,. is the detection significance 
in the x-y energy band). Sixty-one sources were detected only 
in the 0.5-8.0 keV band, while 176 are exclusive to the 0.9-8.0 
keV band and would have otherwise remained hidden because 
of the diffuse emission. This procedure led to a total number of 
1035 X-ray sources in the entire field, which we analyze here. 

3.2. Photon extraction 

Even with the high spatial resolution of the Chandra AClS-1 
camera, the high source density in Trumpler 16 , source photon 
extraction is not an easy task. Although circular regions would 
contain a relatively large fraction of the PSF for almost all 
source photons, the extended wings of the PSF mean that very 
large regions would be needed, incuiTing in the risk of contam- 
ination from nearby sources. Moreover, the resulting inclusion 



of a large number of background events would reduce the sig- 
nal to noise of weak sources. On the other hand, extraction from 
regions that are too small may reduce the photon statistics for 
further spectral and timing analysis. To add ress these issues, w e 
decided to use ACIS Extract (AE) v3.79 (JBrooset al.ll2002h . 
an IDL based package that makes use of TARAI'H CIAO and 
FTOOL:^ software. 

This task reduces the problem of accounting for non- 
Gaussian shapes of the local PSF by calculating the shape of 
the PSF model at each individual source's position. For some 
sources the background level is affected by the extended PSF 
wings of the bright sources (77 Carinae and WR25) in the FOV. 
AE computes source background locally, by defining back- 
ground extraction regions as circular annuli with inner radii 1 . 1 
times the maximum distance between the source and the 99% 
PSF contour, and outer radii defined so that the regions con- 
tain more than 100 "background" events. In order to exclude 
contamination of the regions by nearby sources, background 
events are taken from a "Swiss cheese" image that excludes 
events within the inner annuli radii of all the 1035 sources. 

AE source extraction was performed using a PSF model 
that contains a specified fraction of source events (fpsp)- 
Generally, we choose fpsF=90%, and computed the contours 
from the PSF for a mono-energetic source with E = 1 .49 
keV For 9.9% of the sources in the denser parts of the 
Trumpler 16 field this fraction was reduced to avoid contam- 
ination with other nearby sources, in the most extreme cases 
down to fpsF~50% (just 3 sources). 

Following AE science hints, we then refine the initial 
source positions computed by PWDetecO by correlating the 
source images with the model of local PSF computed by AE 
libraries. This procedure was only used for those sources ly- 
ing at off-axis larger than 5 arcmin (316 sources), while for 
the rest of the source (719 sources) we simply adopt mean 
photon positionO AE also estimates local background spec- 
tra, computes redistribution matrix files (RMFs) and auxil- 
iary response files (ARFs), constructs light curves, performs 
Kolmogorov-Smirnov variability tests, and computes photom- 
etry in 14 different energy bands. Results of AE procedure 
appears in Tabl^U which lists source number in column (1), 
name according to CXC naming convention^ (2), sky po- 
sition (R.A. and Dec. J2000) (3,4) with relative uncertainty 
(5), off-axis angle (0) (6), significance of the detection (Sig.) 
from PWDetect analysis (7), the source extraction area (8); 
the PSF fraction within the extraction area, assuming E=1.49 
keV (9); the background-corrected extracted source counts in 
the 0.5-8.0 keV band (NetCnts) (10); the count rates (CR 
= NetCnts/Exptime/PSFfiac) in three spectral bands: 0.5-8.0 
keV, 0.5-2.0 keV and 2.0-8.0 keV, (11-13); source photon 
quantiles at 25, 50 and 75% percent in columns 14 to 16 
(see sub-section l3.3b . the median photon energy (Ex) in (17). 



** http://asc.harvard.edu/ciao/download/doc/expmap_intro.ps 
' See reasons of this choice in Albacete Colombo et al. ( 2007)). 
'" Afterglow is defined as the residual charge from the interaction of 
a cosmic ray with the CCD. If afterglow events are not removed from 
the data, they can result in the spurious "detection" of faint sources. 



' ' http://www.astro.psu.edu/xray/docs/TARA/ 
'' http://heasarc.gsfc.nasa.gov/docs/software/ftools/ 
'^ PWDetect assumes a symmetric PSF. 

'* Please follow Acis-Extract technical procedures at website 
http://www.astro.psu.edu/xray/docs/TARA/a e_users_guidei 
'^ http://cxc.harvard.edu/cdo/scipubs.html 



Table 1. Trumpler 16 X-ray source catalog (see the electronic version for the complete table). 
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1.01 


-1.05 


-1 





8 


104345.36-593948.7 


10:43:45.37 


-59:39:48.78 


0.36 


8.89 


7.30 


731 


0.89 


54 


0.699 


0.498 


0.201 


0.11 


0.10 


0.31 


1.25 


-0.34 


-2 


0) 


9 


104345.37-593847.5 


10:43:45.37 


-59:38:47.54 


0.40 


9.40 


5.83 


970 


0.91 


45 


0.580 


0.227 


0.353 


0.14 


0.11 


0.49 


1.36 


-2.29 


-1 


&: 


10 


104345.44-594158.9 


10:43:45.44 


-59:41:58.95 


0.44 


8.11 


5.31 


564 


0.90 


25 


0.321 


0.150 


0.171 


0.14 


0.15 


0.43 


1.63 


-0.11 


-1 


> 


11 


104346.39-594929.8 


10:43:46.39 


-59:49:29.82 


0.37 


9.55 


4.72 


474 


0.74 


41 


0.634 


0.550 


0.084 


0.11 


0.10 


0.35 


1.22 


-0.74 


-1 


>< 


12 


104348.15-594924.4 


10:43:48.15 


-59:49:24.41 


0.29 


9.32 


12.27 


927 


0.91 


132 


1.677 


1.276 


0.402 


0.12 


0.10 


0.35 


1.28 


-1.24 


-2 


i-S 

CO 


13 


104349.40-594456.3 


10:43:49.41 


-59:44:56.36 


0.28 


7.45 


13.68 


418 


0.90 


89 


1.138 


1.118 


0.020 


0.06 


0.06 


0.56 


0.98 


-0.96 


-3 




14 


104350.13-594552.6 


10:43:50.14 


-59:45:52.69 


0.47 


7.54 


5.11 


427 


0.90 


15 


0.203 


0.248 


0.044 


0.08 


0.07 


— 


1.01 


-0.36 


-1 


% 


15 


104350.71-593744.4 


10:43:50.71 


-59:37:44.45 


0.42 


9.49 


5.21 


905 


0.89 


41 


0.539 


0.414 


0.125 


0.12 


0.11 


0.47 


1.33 


-0.78 


-4 


CD 


16 


104350.89-595031.2 


10:43:50.90 


-59:50:31.25 


0.39 


9.76 


5.26 


1073 


0.90 


23 


0.302 


0.164 


0.138 


0.16 


0.09 


0.68 


1.15 


-2.67 


-4 


a. 


17 


104351.09-594024.6 


10:43:51.10 


-59:40:24.67 


0.38 


7.97 


5.42 


479 


0.89 


28 


0.361 


0.239 


0.122 


0.11 


0.09 


0.36 


1.20 


-0.31 


-1 


% 


18 


104351.64-594525.1 


10:43:51.64 


-59:45:25.10 


0.40 


7.25 


7.05 


376 


0.90 


31 


0.400 


0.253 


0.147 


0.14 


0.16 


0.39 


1.72 


-3.15 


-2 


B 


19 


104351.87-594035.6 


10:43:51.88 


-59:40:35.60 


0.43 


7.81 


5.87 


443 


0.89 


7 


0.102 


0.102 


0.001 


0.13 


0.09 


— 


1.21 


-0.34 


-2 




20 


104352.13-594802.0 


10:43:52.13 


-59:48: 2.00 


0.40 


8.15 


5.23 


528 


0.90 


30 


0.388 


0.239 


0.149 


0.12 


0.09 


0.37 


1.16 


-0.46 


-4 


Vi 

% 


21 


104352.25-594157.6 


10:43:52.25 


-59:41:57.68 


0.39 


7.28 


5.31 


357 


0.89 


25 


0.325 


0.301 


0.023 


0.07 


0.08 


0.93 


1.07 


-0.16 


-4 


22 


104352.48-593920.9 


10:43:52.48 


-59:39:20.96 


0.23 


8.35 


19.14 


592 


0.90 


176 


2.250 


1.401 


0.849 


0.13 


0.16 


0.39 


1.70 


-0.50 


-2 


3' 

H 


23 


104354.13-594145.2 


10:43:54.14 


-59:41:45.24 


0.37 


7.11 


4.81 


357 


0.90 


25 


0.325 


0.101 


0.224 


0.09 


0.15 


0.38 


1.63 


-0.17 


-1 


B 
'a 


24 


104354.20-593805.2 


10:43:54.21 


-59:38: 5.24 


0.41 


8.93 


5.21 


737 


0.89 


22 


0.290 


0.196 


0.094 


0.15 


0.12 


0.72 


1.43 


-0.52 


-1 


25 


104355.11-593624.2 


10:43:55.11 


-59:36:24.23 


0.30 


10.04 


13.22 


1205 


0.90 


177 


2.246 


1.378 


0.869 


0.13 


0.15 


0.37 


1.59 


-4.00 


-2 


« 


26 


104355.14-594750.4 


10:43:55.14 


-59:47:50.42 


0.35 


7.72 


7.89 


451 


0.90 


44 


0.567 


0.470 


0.097 


0.07 


0.08 


0.50 


1.07 


-1.07 


-3 


S 


27 


104355.47-594253.5 


10:43:55.48 


-59:42:53.59 


0.29 


6.69 


13.04 


278 


0.89 


51 


0.664 


0.357 


0.306 


0.14 


0.17 


0.53 


1.74 


-3.51 


-2 


g. 


28 


104355.56-594923.0 


10:43:55.57 


-59:49:23.07 


0.30 


8.57 


8.43 


629 


0.90 


77 


0.989 


0.807 


0.182 


0.09 


0.08 


0.35 


1.13 


-0.50 


-2 


n 


29 


104356.25-594936.4 


10:43:56.25 


-59:49:36.48 


0.33 


8.65 


5.36 


648 


0.89 


36 


0.471 


0.335 


0.136 


0.14 


0.09 


0.39 


1.21 


-0.41 


-1 




30 


104356.82-594236.0 


10:43:56.82 


-59:42:36.03 


0.46 


6.57 


4.90 


257 


0.89 


12 


0.158 


0.139 


0.019 


0.12 


0.12 


— 


1.41 


-0.41 


-1 


_ s 



Notes: Column labeled with refers to the off-axis angle measured in arcmin from the aim point of the observation. Sig. is the significance of the source in number of sigma over background. 
Flag Id = -1: sources detected only in the 0.9-8.0 keV band, flag Id = -2 and -3 correspond to sources detected in both bands, but Sigo.9-8.0 > Sigo.5-8.0 and Sigog-so < Sigos-so. respectively. Flag 
Id = -4 refers to sources detected only in the 0.5-8.0 keV band. 
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Column 18 is the log(P|js) Kolmogorov-Smirnov probability 
of non-variability (see section §|5]l, and in (19) are flags from 
PWDetect detection code. 



3.3. X-ray hardness ratios 



A commonly used tool to explore the spectral properties of 
sources with low phot on statistics is the hardness ratio (e.g. 
Schulz et al.L !l989; Pre stwich et al.Ll2003b . In this conventional 
method, the full energy range is divided into two or three sub- 
bands and the detected source photons are counted separately 
in each band. Most popular definitions for a single hardness 
ratio (HR) exists on basis of only two energy sub-bands: /- 
HR=H/S or a- HR=(H-S)/(H+S). By these definitions, HR is 
very sensitive to small changes (i.e. statistical fluctuations) in 
the number of photons falling in each band. The requirement of 
total counts in the full energy band is at least 40 photons (just 
36% of our sources satisfy this constraint). Above this limit 
HR becomes a "reUable method" to estimate the real hardness 



of sources (lAlbacete Colombo et al.L l2007h 



An improved met hod to resolve this limitation is based on 
the quantile analysis ( Hong et al.l 120041) . Instead of working 
with predetermined energy bands, we determine the energy E^ 
below which the net counts is x% of the total counts of the 
source. We define quantile Q^ as: ^ ~_^"" , where in our study 
Emin=0.5 keV and Eniax=8.0 keV. We compute median Qso val- 
ues and quartiles Q25 and Q75 and give values in Table[T] A 
minor inconvenience of this method is that for a given spec- 
trum, various quantiles cannot be considered independent vari- 
ables, unlike the counts in different energy bands. However, 
Hong et al.l (120041) overcome this problem by considering the 
log(Q5o/(l-Q5o)) vs. 3(Q25/Q75) plane. Based on an extensive 
set of simulated spectra, they predict loci of models in this 
plane. We use a set of absorbed thermal models with plasma 
temperatures of 0.2, 0.5, 1, 2, 4, 10 keV and Nh equal to 10^°, 
10^1, 0.4x10^2, 10^2, 4x10^2 and 10^^ cm^^. Note that the spec- 
trum changes from soft to hard as one goes from left to right in 
the diagram (see Fig|2]i. 

In a statistical sense, a comparison between sources in 
the grid models suggest typical Nh and kT values distributed 
around O.SxlO^^ cm ^ and 1.5 keV, but dispersed within grid 
models of Nh=0.4-1.0x1022 cm^^ ^^^j kT=l-3 keV, respec- 
tively. The absence of a systematic difference between the po- 
sition of sources with and without a 2MASS identification 
on the hardness-ratio plane is likely to a combination of two 
observational bias: / related to the limiting magnitude of the 
2MASS photometry (Kj,'™' -14.3 mag.); // the X-ray sensitiv- 
ity of the Chandra data (f|('™' -lO"'"*'^ erg s"^ cm"^) However 
is not clear, from the observational point of view, i.e. throgouh 
the near-IR and X-ray flux source distributions, that a fraction 
of faint stellar X-rays sources would not have 2MASS coun- 
terparts. Deeper near-IR and X-ray observations are needed to 
unveil this issue. 



0.01 

_kT (keV) = 0.2 




-1.0 -0,5 

log,o(Q5o/(1-Q5o)) 



Fig. 2. Quantile Color-Color Diagram (QCCD). The energy 
scale in the top X-axis shows the median energy values (Qso). 
The grid patt ern represents the location of simulated spectra in 
the diagram (iHong et al.l l2004l) . Note: Filled circles and open 
boxes refer to X-ray sources with and without 2MASS coun- 
terparts (see next section), respectively. Small dots (red) cor- 
respond to sources with unreliable quantile values affected by 
poor photon statistics (< 10 photons). 



4. Optical and near-IR analysis 

4.1. Counterparts 

The large amounts of gas, dust and selective extinction of the 
region, combined with the absence of available deep (V>19 
mag) optical photometry, prevent finding optical counterparts 
for most of our X-ray sources. We consider it appropriate only 
to use opti cal data to identify the mas sive star population (typi- 
cal V <12. iMassev & Johnsoiild 1993b ) of the region (i.e. Wolf- 
Ray et, O- and early B-types). We use a recent c ompendium o f 
massive stars in the Carina Nebula presented by SmithI (l2006h . 
Out of a total of 60 stars, 44 lie within the 0.0823 deg^ FOV of 
our X- ray observation. We a dded the binary (05.5V + 09.5V) 
F015 (iNiemela et all |2006'), not included in the list of lSmith 
(l2006l) . Our final list of 45 hot massive stars is comprised of: 1 
LBV, 1 Wolf-Rayet, 21 O-type stars and 22 early B-type stars. 
With a cross-identification radius of 3 arcsec, a total of 28 X- 
ray sources were identified. All the O-type stars were identified 
as X-ray emitters, while just 7/22 (~31%) of early B-type (SpT. 
between BOV to B1.5V) have detectable X-ray emission. The 
list of massive stars with X-ray counterparts and the discussion 
of their X-ray properties is presented in §|8] 

Apart from the massive stars, our X-ray population is essen- 
tially composed of low-mass stars. In such a young cluster like 
Tr 16, low-mass stars are often optically invisible, being em- 
bedded and/or obscured by high amounts of gas and dust. We 
partially solve this problem by use of the near-IR observations, 
on which the impact of dust extinction is reduced. We adopt 
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J (1.25 jum), H (1.65 fim) and Kj (2.17 yum) photometry from 
the Two Micron All Sky Survey (2MASS) Point Source Catalog 
(PSCH 2MASS is complete to magnitudes of 15.8, 15.1 and 
14.3 mag in the J, H and K, bands, respectively. We restricted 
our photometry analysis to 2MASS sources with flag quality 
A, B, C or D in at least one of the three magnitudes (see ex- 
planation in the 2MASS Afl-Sky Data Release User's Guide). 
With this restriction 20 sources were removed from our initial 
list of 5938 sources in the ACIS FOV of our observation. This 
leaves a total of 5918 2MASS object^ 

We began by cross-identifying our X-ray source list with 
the 2MASS catalog. Identification radii. Rid, were chosen to 
limit the number of spurious identifications due to chance 
aUgnments, Nchance, and at the same time to include a large 
number of the true physical associations, Ntme. Identification 
radii used were 1.0, 1.5,2.1, and 2.7 arcsec adopted for the four 
different off-axis angle ranges: [0-2), [2-4), [4-7) and > 7 ar- 
cmin, respectivel}!^. Results of the final identification are pre- 
sented in § 14.31 and shown in Table|2l The first seven columns 
are: X-ray source number; 2MASS nomenclature for identified 
sources; the ofi'set between the two positions; J, H, K, magni- 
tudes; photometry quality flag (Ph.Q); confusion flag (Cont). 
A total of 660 X-ray sources out of the 1035 in our list were 
identified with 665 2MASS objects. Five X-ray sources (#96, 
#382, #401, #816 and #1034) were identified with two 2MASS 
counterparts each. After a visual inspection, we kept only the 
closer counterparts. Thus the final list of near-IR counterparts 
consists of 660 entries. 

We estimate the expected n umber of extragalactic sources 
in our detection list by following lFlaccomio et al.l (|2006 j proce- 
dure. We consider the ACIS count-rates of non-stellar sources 



in the Chandra Deep Fie ld North (CDFN, [Alexander et al. 



2003t iBarger et al.L l2003l) and estimate absorption corrected 



count-rates assuming Nh = 5x10~' cm"- (from Av~3.6, see 
§ 14.31 ) using P IMMS and assuming p ower-law spectra with in- 
dex 1 and 2 ( IGiacconi et al.l l200lh . We then compare these 



count rates with upper limits taken at random positions in the 
ACIS FOV. For F between 1 and 2 we obtain 72 to 95 ex- 
pected extragalactic sources. Given the intrinsic near-IR fluxes 
of these sources and the absorption toward Trumpler 16 , they 
are expected to be among the 385 without NIR counterparts (cf. 
Flaccomio et al. 2006). This means that no more than 18 to 24% 
of the unidentified X-ray population is related to extragalactic 
sources. 

4.2. Unidentified X-ray sources 

A large population of young stars, proto-stars, deeply embed- 
ded in dense cir cumstellar gas and dust should be present in the 
Carina Nebula (ISmith et al.Ll2003h . However optical and near- 
IR counterparts of YSOs are difficult to detect. Fortunately, X- 
ray emission is expected in young stellar objects along all their 



'* See http://www.ipac.caltech.edu/2mass' 

" We note that the 2MASS catalog appears to have a "hole" around 
TyCarinae (~ 1.5 arcmin radius). 

"* The ado pted radii were computed follow ing technical procedures 
presented in lAlbacete Colombo et alj ( l2007h 



initial phases (IMontmerle et al.L 120001) . X-ray photons easily 
escape from dense circumstellar material, where absorption 
process becomes important mainly for energies below 1.2 keV 
(Morrison & McCammon, 1983b). This makes hard X-ray en- 
ergies the most appropriate "window" to detect counterparts of 
deeply embedded young sources. 

Of the 375 X-rays sources without near-IR counterparts, 
just a small fraction (~20%) is expected to be extra-galactic 
contamination (see 94.11 1. We are thus dealing with about ~300 
candidate young (first stage) low-mass stars, highly obscured 
by circumstellar material. In our data these sources typically 
have lower X-ray photon statistics than those with near-IR 
counterpart, i.e. ~8 vs. 28 average photons, respectively. In 
X-rays, no quantitative differences in the median energy and 
spectral quantiles were found for the X-ray sources with and 
without near-IR counterparts. 

4.3. Near-IR properties of identified X-ray sources 

We now investigate the near-IR properties of the X-ray 
sources. For this purpose we restrict our analysis to sources 
with high quality photometry (Ph.Q=AAA) and no confusion 
(Co«f.=000). With these requirements the total number of IR 
sources in the ACIS FOV is reduced from 5918 to 2178. We 
have also set a further requirement on near-IR counterparts of 
X-ray sources, that their J, H, and Ks magnitude eiTors be all 
<0. 1 mag. All these requirements yield 367 X-ray sources with 
good near-IR counterparts, out of the original 660. 

Figure[3] shows the J-H vs. H-Ks color-color (CC) dia- 
gram for these AAA-flagged sources. We a lso plot for com- 
parison the MS (IKenyon & H artmann. 1995), the Classical T- 
Tauri Stars (CTTS) locus of [Mever et al.. (.1997.) . and redden- 
ing vectors starting from these loci and with slope (Aks/E(H - 
Ks)=0.I 2 5) co rresponding to the extinction law given by 
HansonI (l2003h . Trumpler 16 members with purely photo- 
spheric emission should he in this reddening band. Otherwise, 
Young stellar objects (YSOs), such as Classical TTauri and 
Herbig Ae/Be stars, because of the NIR excess emission orig- 
inating in the inner parts of their circumsteUar disks, are often 
found to the right of this band, i.e. in the CTTS locus. Fifty-one 
(out of 367) X-ray sources, i.e. likely Trumpler 16 members, 
have colors consistent with the (reddened) CTTS locus. This 
means a fraction of 51/339 (28 OB stars were discarded) ~ 15% 
of all (low-mass) identified X-ray sources in the CC diagram. 
Of all 51 sources with intrinsic K-excess i.e. disk-star systems, 
eleven (Src-Id: 36, 41, 209, 230, 773, 966, 993, 996, 1002, 1003 
and 1009) appears below the CTTS vectors, but with intrinsic 
bright Ks magnitudes, as is shown in both panles of FiglH They 
are probably intermediate- to high-mass young stars with an 
intrinsic K-band excess that would be produced by massive ac- 
cretion disks and/or extended envelopes surrounding massive 
YSOs. If confirmed, they will contribute for about 40% of the 
total massive star population of the Trumpler 16 region. 

FigureHl-left shows the K, vs. J-K, color magnitude (CM) 
diagram for the same stars plotted in Fig. [3] We also show 
for reference the expected cluster locus: the intrinsic Ks mag- 
nitudes and J-Ks colors for stars earlier than B5V were 
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Table 2. Near-IR counterparts of Trumpler 16 X-ray sources. The complete version is available in the electronic edition of A&A. 



Nx 


2MASSJ-I- 


Offset (") 


J mag. 


H mag. 


Ks mag. 


Ph.Q 


Cent. 


Av 


Mass Notes 


1 

2 
3 


10433859-5938306 


2.33 


14.57+0.03 


13.21+0.04 


12.63+0.03 


AAA 


000 


7.93 


2.05 


10434145-5942245 


0.30 


14.18+0.03 


13.20+0.02 


12.82+0.02 


AAA 


000 


4.59 


2.24 


4 


10434126-5941002 


2.62 


15.47+— 


15.31+0.14 


14.60+0.14 


UBB 


000 


— 


— 


5 
6 

7 


10434295-5944080 


1.77 


15.63±0.08 


14.31+0.06 


13.27+0.04 


AAA 


000 


15.27 


1.25 K-excess 


10434401-5948177 


0.76 


8.70±0.02 


8.51+0.04 


8.47+0.02 


AAA 


000 


1.38 


19.2 09.5V 


8 


10434538-5939468 


1.92 


15.90+-- 


14.88+— 


14.90+0.16 


UUC 


000 


— 


— 


9 


10434536-5938471 


0.36 


16.02+0.10 


14.04+0.04 


13.02+0.04 


AAA 


000 


14.90 


1.01 


10 


10434525-5941567 


2.61 


14.37±0.06 


12.99+0.05 


12.18+0.04 


AAA 


000 


12.05 


2.16 K-excess 


11 


10434659-5949292 


1.68 


13.69±0.05 


12.68+0.05 


12.25+0.04 


AEA 


cOc 


— 


2.38 


12 


10434809-5949246 


0.50 


13.51±— 


13.09+0.07 


12.85+0.05 


UAA 


Occ 


— 


— 


13 


10434937-5944549 


1.45 


12.67+0.02 


12.20+0.03 


12.06+0.03 


AAA 


000 


0.82 


4.89 


14 


10435007-5945530 


0.59 


15.56+0.05 


14.59+0.03 


14.19+0.07 


AAA 


000 


4.31 


1.30 


15 


10435085-5937437 


1.28 


14.86+0.05 


13.74+0.05 


13.40+0.05 


AAA 


000 


3.76 


1.86 


16 


10435088-5950307 


0.47 


12.20+0.02 


11.99+0.03 


11.88+0.02 


AAA 


000 


0.28 


6.22 


17 


10435123-5940243 


1.06 


14.60+0.05 


13.38+0.04 


12.88+0.04 


AAA 


000 


6.57 


2.03 


18 


10435132-5945239 


2.64 


15.54+0.08 


14.06+— 


13.55+— 


AUU 


cpp 


— 


1.31 


19 


10435191-5940353 


0.40 


16.49+0.16 


15.15+0.09 


14.63+0.11 


CAA 


000 


— 


0.70 


20 


10435186-5948017 


2.02 


14.49+0.05 


13.46+0.06 


12.95+0.05 


AAA 


000 


6.60 


2.09 


21 


10435223-5941574 


0.29 


14.83+0.05 


14.46+0.08 


14.35+0.11 


AAA 


000 


— 


1.88 


22 


10435230-5939222 


1.87 


13.03+0.04 


11.90+0.06 


11.20+0.04 


AEE 


000 


— 


3.55 Mass-deg. 


23 


10435408-5941463 


1.19 


15.40+0.07 


14.24+0.05 


13.79+0.05 


AAA 


000 


5.47 


1.40 


24 


10435419-5938073 


2.13 


13.03+— 


13.70+0.06 


13.53+0.07 


UAA 


Occ 


— 


— 


25 


10435501-5936242 


0.75 


11.64+0.02 


11.48+0.03 


11.37+0.03 


AAA 


000 


0.42 


8.28 


26 


10435505-5947505 


0.70 


14.18+0.03 


13.81+0.04 


13.70+0.06 


AAA 


000 


— 


2.24 


27 


10435545-5942531 


0.50 


14.93+0.05 


14.11+0.05 


13.84+0.06 


AAA 


000 


2.37 


1.81 


28 


10435557-5949226 


0.39 


13.44+0.03 


12.65+0.03 


12.24+0.03 


AAA 


000 


5.08 


2.45 


29 


10435606-5949351 


1.97 


14.53+0.05 


13.83+0.05 


13.60+0.06 


AAA 


000 


1.87 


2.07 


30 


10435684-5942364 


0.45 


14.70+0.05 


13.57+0.03 


12.88+0.04 


AAA 


000 


9.63 


1.97 K-excess 



Column 3 ("Offset") is the offset between X-ray and near-IR counterpart. Ph.Q refers to the 2MASS photometric quality flags for the J, H 
and, Kj bands: "A" to "D" indicate decreasing quality of the measurements, "U" that the value is an upper limit. The next column refers to the 
contamination and confusion flag: For further analysis we considered only sources unaffected by known artifacts, i.e. Cont. = 000 (see 2MASS 
documentation for details). Masses are given in solar units, and the last column contains information presented in §|4] Note: the "Mass-deg." 
flag indicates mass degeneracy according to a Siess-based Jmag-Mass calibration. 



taken from the MS c alibration of iKnodlsedeii (120001) and 
Bessell & Bretli (Il989h . For later spectral types (masses be- 
tvyeen 0.1 a nd 7 M p), we adopted the 3 Myr isochrone from 
Siess et al.l (120001) . convert ed to the observation a l plan e us- 
ing the cahbration given bv lKenvon & HartmannI (1 1995b . The 
adopted MS and 3 Myr isochrone overlap satisfactorily. 

In order to estimate the typical visual absorption of clus- 
ter members, we computed the distance of each X-ray source 
from the cluster locus along the reddening direction on the 
KsVsJ-H plane (see Fig.|4]-right). Resulting Ay values for in- 
dividual sourceo are listed in column 9 of Table|2] Note that 
for 11.2<Ks <12.8 the absorption cannot be constrained be- 
cause the reddening vector intersects the cluster locus more 
than once. In both CM diagrams, some X-ray sources, of the 
order of 20, lie to the left of or relatively close to the unred- 
dened cluster loci. These are likely to be foreground MS stars 



" Individual photometric errors of up to 0. 1 mag at Ks and J-H color 
errors of ~0.14, could results for Ay errors up to ~0.7 mag. 



and thus their Av values either cannot be computed or are close 
to zero. In a statistical sense, the median Ay value of OB stars 
(Ay =2.0+0.8 mag.) is lower than that computed for low-mass 
stars (Ay =3.6±2.4 mag., considered to be the typical absorp- 
tion of the cluster). Obviously, the above estimates depend 
on the reliability of the assumed cluster locus and on the as- 
sumption that disk-induced excesses do not significantly af- 
fect the J and H magnitudes. To caution of possible contam- 
ination and/or anomalous J magnitude, but in particular for 
the H-band, the median Ay for the low-mass star population 
was estimated by discarding sources with intrinsic Ks excesses 
(i.e. labeled with asterisks in figures [3] and |4|. We must note 
that: /- the dispersion along the J-Ks axis (see FigHl-left) in- 
dicates differential absorption of the region, and translates in a 
broad Ay distribution, with a spread Icr ~2.4 mag. //- the dif- 
ference between median Ay values of low-mass and OB stars 
suggests a clearing effect of strong winds and radiation field of 
massive stars on their suiTounding environment. This conclu- 
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Fig. 4. CM diagrams of the Trumpler 16region. Symbols as in figurelH The two parallel curves indicate the expected cluster loci 
for the assumed distance and no reddening, and for the mean reddening Ay =3.6 mag. Masses are indicated on the left side of 
the unreddened isochrone. Identified OB stars show Kj magnitudes below 11. Stars with Kj excess show typically J-Kj > 1. 
Note: The limitation imposed by the 2MASS photometry (K, ~14.5 mag) corresponds to a mass limit of 1.1 M©. Labels near 
OB stars refer to Simbad names (left) and respective Spectral types (right). The fact that stars Trl6 - 244 (04If) and HD93 129AB 
(WN6ha-i-04f) appear above the 120 M© track does not mean that they are more massive, since their evolutionary stages cannot 
be compared with a MS isochrone. //Carinae was discarded because of the bad 2MASS photometry (Phot.qual.='DDD', see 
text). Sources labeled with numbers indicate stars (with masses above 12 M©) probably surrounded by disks. 



sion was also reached by JAlbacete Colombo et alj (12007!) for 
massive stars of the Cyg OB2 region. However, compared to 
the Cyg OB2 region, the fraction of disk-stars members in the 
Trumpler 16 region is about four times larger. 

Finally, we use the 2MASS J-band magnitudes to obtain an 
estimate of stellar masses for 510 counterparts (of a total 660 
identified stars) with J-band Phot.Qual.-' N to 'D'. We com- 
pute the mass vs. J mag relationship appropriate for the cluster 
mean age (3 Myr), distance (DM= 1 1 .78 mag) and extinction 
(Av =3.6 mag), this latter obtained as describ ed above for the 
cluster locus in the CM diagramo- We use Siess et al.l ( 12000) 



models to compute PMS tracks of low- and intermediate-mass 
stars (masses < 7 Mq). We interpolated the J mag - mass rela- 
tion, using the J magnitude vs. mass relationship at 3 Myr (see 
column 10 of Table|2]l, to compute individual masses of staro- 
Unfortunately, the relation mass-J mag degenerates in the mass 
ranges 0.2-0.65 M© (3 sources) and 2.7-4.52 M© (21 sources). 
We indicate these sources with "Mass-deg." flag in column 1 1 
of Table|2] and give mean mass values for sources lying in these 
two ranges, i.e. 0.42 and 3.55, respectively. Computed masses 



-" The choice of the J band is justified because (i) in the presence of 
disk excesses the J-band is the most representative of the photospheric 
emission i.e. least affected and (ii) the mass ranges in which the mass- 
luminosity relationship is degenerate are narrower than for a similar 
relationship in the H and Ks bands. 

-' This method suffers of photometry inaccuracy, distance and age 
spread. Mass values should be adopted carefully and should not be 
considered to measure the slope of the mass function. 



over 7 M© are potentially affected by large uncertainties in the 
extrapolated J mag - mass relation, and these values were ex- 
cluded from further analysis. In Table |2] we give masses for a 
total of 510 stars, of which 410 range between 0.65 and 2.52 

Mo. 

5. X-ray variability 

PMS stars have high levels of X-ray activity that are com- 
monly attributed to a "scaled up" solar-like corona formed 
by active regions. X-ray variability over a wide range of 
time scales is common in all magnetically act i ve stars (e.g . 



Ipgig elson & Montmerleill999l:lFavata & Micelall200llGudelL 
2004). On long time scales, this includes rotational modulation 
of active re gions, their emergence and evolution and m agnetic 



cycles (e.g. iMarino et all 120031 : iFlaccomio et all 12005.) . Most 
of the observed variations have short time-scales (~ hours), 
however, and can be attributed to small scale flares triggered 
by magnetic reconnection events. 

We first investigated X-ray variability in our sources us- 
ing the non-binned one-sample Kolmogorov-Smirnov (KS) test 



(IPress et al.Lll992h . This test compares the distribution of pho- 



ton arrival times with that expected for a constant source. The 
test was applied to photons in the source extraction regions, 
which also contain background photons. Given that the back- 
ground was found to be constant with time (§2.1), the results, 
i.e. the confidence with which we can reject the hypothesis that 
the flux was constant during our observation, can be attributed 
to the source photons. TablefT] column 18, reports the logarithm 
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Fig. 5. Light curves (in the 0.5-8.0 keV band) showing different variability scenarios occurring among our 77 variable sources. 
Bin size is 900 sec. The source number is indicated in each panel. Upper panels show a flare-like behavior with very quick rises 
and decay phases of hours, while bottom panels are slowly modulated variability. A probable flare-i-modulation variability is 
occurring in the last panel. 



of the KS-test significance with values < -4 truncated at that 
value: sources with log(PKs) <-3.0 can be considered almost 
definitely variable as we expect at most one of the 1035 sources 
(i.e. <0.1%) to be erroneously classified as variable. Seventy- 
seven X-ray sources (xlA% of the total) fall in this category. 
Fifty-five sources with -2.0 < log(PKs) < -3.0 can be considered 
as likely variable, although about half of them suffers of low 
photon statistics. These numbers of sources are lower limits 
to the total number of variable sources in the region for sev- 
eral reasons: i) most of the observed variability is in the form 
of flares, i.e. events that are shorter than our observation and 
with a duty-cycle that may be considerably longer ( Wolk et al., 
2005h : ii) the sensitivity of statistical tests to time variability of 
a given relative amplitude depends critically on photon statis- 
tics (see Albacete Colombo et al. (2007)). Hereafter, we con- 
sider variable those 77 sources with log(Pks)<-3 . 

To get a more accurate description of the detected vari- 
ability, we extracted binned light-curves for each of the 77 
variable sources in the region. We adopt a bin length of 900 
seconds, a compromise between bins that are long enough to 
reach a good signal-to-noise ratio per bin for most sources and 
sufficiently short to resolve the decay phase of typical flares. 
Since the background of our observation is both low (negli- 
gible for many sources) and constant in time, we did not ap- 
ply any background subtraction to the presented light-curves. 



In Fig |5] we show examples of the dififerent behavior among 
the light-curves of variable sources. Source #254 like others 
(#35, #270, #438, #468, #480, #596, #761 and #975) experi- 
ence "impulsive" flares with very quick rises and decay phases 
of only a few hours. Others (sources #503 (#36, #41, #136, 
#164, #260, #489, #523, #731 and #811) show longer (2 to 10 
hours) flares. In several instances a second impulsive event is 
visible during the exponential decay of a previous flare (e.g. 
sources #36, #52, #87, #196, #564, #600, #620 and #623). The 
case of source #696 is a combination of both variability types, 
with two consecutive flares. Other sources like #714 and #651 
(e.g. #27, #251, #271, #793, #839, #890, #904 and #980) show 
light-curves that, instead of showing typical flares, are charac- 
terized by slow continuous rises or decays that might be ex- 
plained by rotat ional modulation of no n-homogeneously dis- 
tributed plasma (IFlaccomio et al.l l2005b . Finally, light-curves, 
like those observed for sources #524, #71 and #933, seem to 
be related to a combination of flare like activity and rotational 
modulation. 

5.1. Variability in massive stars 

Because X-ray emission from O stars, which is believed to 
be unrelated to solar-like magnetic activity, comes from the 
integrated emission from a large number of small shocks 
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Fig. 3. JHKj color-color diagram from high quality 2MASS 
photometry. Filled circles and dots refer to X-ray detected and 
undetected 2MASS sources, respectively. Asterisks show the 
source population with intrinsic Kj excess. Identified 2MASS 
source known as O- and early B-type stars are indicated by 
SQUARES. The main seq uence is shown for reference. We also 
show the CTTS locus of lMever et al.l (119971) and reddening vec- 
tors {dashed lines) with length corresponding to Ay =13 mag. 
Note: i-: the 2MASS photometry of zyCarinae is severely af- 
fected, with Phot.Qual.-'DDD' ; ii-: the peculiar position of 
the red box with H-K, <-0. 1 corresponds to the massive binary 
CPD-592628AB (SpT; O9.5V+B0.3V). 



randomly occurring in their st rong winds (IFeldmeier et al. 



1997t lOwocki & Cohenl 119991) . on average, global X-ray 
variabiUty is not expected to occur. However, it is surpris- 
ing that three (out of 28) massive stars, namely Tr 16-11 
(B1.5 V, source#136), Trl6-5 (Bl V, source #489) and the bi- 
nary HD 93205 (03.5V((f))H-08V, source #242), are signifi- 
cantly variable, with log(PKs) values lower than -3. The origin 
of the observed flare-like variability in the first two sources (B- 
type stars) is probably coronal activity of unresolved late-type 
compan ions. This hypothe sis was also proposed for Trl6-ll 
itself bv lEvans et al.l(l2003h . 

Fig|6] shows the X-ray light-curve of the HD 93205 binary 
system. Two different processes may be acting simultaneously 
to explain the observed variability: 

i- The decrement of count r ate agrees with phas e-locked X- 
ray variability reported by lMorrell et al.l (1200 ih . The X-ray 
count rate decreases from ~ 0.055 to 0.033 cts/s (about 
35%) in about 88 ksec (~1 day) of continuous observa- 
tion, i.e. about 16% of the orbital period. This is in agree- 
ment with X-ray emission from colliding winds contribut- 
ing most of the total detected emission. However, we can- 
not discard that magnetically channeled wind streams, at 
the 03.5V((f)) primary, collide with the cool and dense 
postshock plasma at the magnetic equator (e.g. Orionis-C, 
Gagne et al.l 12005 ah . A detailed spectral and time (phase- 
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Fig. 6. X-ray light-curve of the massive binary 
(03.5V((f))-H08V) HD 93205. This light-curve shows: i- 
a uniform decay probably related to the orbital motion of the 
system, with changing absorption Nh of the colliding-wind 
region along the line of sight, plus ii- a probable short-term 
variability (see text). 

resolved) study should be follow to discern between these 
two possibilities. 
//- The observed short time X-ray v ariability (log(PK s)=- 



3.07). While dynamical instabilities (iLuo et al.Lll99' 



JOF 



in 

the Colliding Wind Region (CWR) seems to be more ap- 
propiate, the magnetic reconnection from an unknown low- 
mass tertiary companion^ cannot be rule out. 

We remark that the origin of X-ray variability in high mass 
stars is out of the scope of this work, it needs for a large X-ray 
data-set. We leave this subject to a further investigation. 



6. Spectral analysis 

In order to characterize the hot plasma responsible for the X- 
ray emission of Trumpler 16 stars, and to estimate their intrin- 
sic X-ray luminosities, we analyzed the ACIS spectra of the 
615 (out of all 1035) sources with more than 20 net photons 
(NET_CNTS), i.e. corrected for local background. Spectral pa- 
rameters for sources with less than 20 net photons are too 
much ill-constrained ( lAlbacete Colombo et al.ll2007b . and thus 
were not determined. Moreover, high local background could 



^- Whenever the wind velocities of two stars are not equal, shocked 
CWR should be subject to the rapid growth of dynamical (Kelvin- 
Helmholtz) instabilities, reaching the radiative cooling state. This limit 
becomes i mportant for massive binaries with typical orbital periods 
<20 days JAntokhin et"aill2004l) 

^^ Variability like HD 93205's has been seen also on theta 20ri- 
A and has been interpreted by Schulz et al. (2006) as evidence of 
binary-induced reconnection. However, Image reconstruction shows 
that hard, short-term variat ions are from an int ermediate-mass tertiary 
companion at ~0.3 arcsec Marvin et alil2002h . 
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affect the reliability of computed spectral parameters. We de- 
fine the fraction fcont as the ratio between BKG_CNTS (num- 
ber of background counts in the source extraction region) and 
NET_CNTS. We accepted spectral fits for sources with fcont<l 
and NET_CNTS>20 photons. Of all 1035 sources, only 563 
satisfy both conditions above, while 119 show fcont<l and 353 
lie in the low statistics regime (i.e. NET_CNTS<20 ph). 

Source and background spectra in the 0.5-8.0 keV band 
were produced with AE (see 33. 2b . along with individual 
"redistribution matrices files" (RMF) and "ancillary response 
files" (ARE). For model fitting, spectra were grouped so to 
have a specified number of events in each energy bin. Grouping 
was tuned to the source statistics and we chose 2, 5, 7, 10, 
and 60 counts per channel for sources with net-counts in the 
following ranges: [20-40], [40-100], [100-200], [200-500], and 
[500-**]. Spectral fitting of back ground-subtrac ted spectra was 
performed with XSPEC vl2.0 dArnaudi |2004 and our own 
shell an d TCL s cripts to automate the process as described in 
IFlacco mio et al.l ([2006). Because background c orrected spec 



Table 3. X-ray spectral fits of Trumpler 16 sources. The com- 
plete version is available in the electronic edition of A&A. 



tra ar e not appropriately handled by C-statistics (IGetman et al 



12005 ^. best-fit parameters for the chosen models were com- 
puted by means of the Chi-Squared (x^) minimization. 

We fit our spectra assuming emission by a thermal plasma, 
in collisional ionization equilibrium, as modeled by the 
APEC code (Smith etal., 2001). Elemental abundances are 
not easily constrained with low-statistics spectra and were 
fixed at Z=0.3 Zq, with solar abundance ratios taken from 



ia 



Anders & Grevessd (119891) . The choice of sub-solar abun- 



dances is s uggested by several X-ray studies of sta r forming re 



gions (e.g. F eigelson et alil2002tlPreibischLl2003h . Absorption 



was accounted for using the WARS model, parameterized by 
the hy drogen column density, Nh dMorrison & McCammoni 
1983al) . In Table|3]we give best-fit parameters (xl, Nh, kT and 



Lx) of the sources. 

Except for 28 massive O- and early B-type stars, we fit 
source spectra with one-temperature (IT) plasma models using 
an automated procedure. In order to reduce the risk of finding a 
local minimum in the x^ spaces, our procedure chooses the best 
fit among several obtained starting from a grid of initial values 
of the model parameters: log(NH) = 21.0, 21.7, 22.0, 22.4, 22.7 
and 23.0 cm^ and kT = 0.5, 0.75, 1.0, 2.0, 5.0 keV. best-fit val- 
ues of log(A^H) < 20.3 cm~^ were truncated at 20.3 for two cases 
(sources #150 and #944) because, in the 0.5-8.0 keV energy 
range, ACIS spectra are insensitive to lower column densities. 
In a similar way, above 10 keV Chandra is not able to discrim- 
inate between such high temperatures. Therefore, 74 best-fit 
values of kT above 8 keV were truncated to that value. They 
are indicated with a flag hard-tae. in Table[3] 

Fig.E] shows the distribution of best-fit log(NH) values for 
the 563 fitted sources. They appear to be normally distributed 
with a median logA^n -21.73 (Nh=5.37x102' cm"^) and a 
FWHM of ~0.4 dex. The log-normal distribution of the sources 
is indicated with the Gaussian curvq^- The computed median 
of Nh (~5. 37 x10^' cm~^) is con verted to a median Ay =3.35 
by use of the IVuong eTail (l2003l) relation: Nh/Av = 1.6x10^' 



N. 


Cnts. 


Stat. 


log(NH) 


kT 


logCU) 


flag 


# 


(ph) 


(xt) 


(cm"2) 


(keV) 


(erg/s) 




1 


30 


— 






30.51 


no-fit 


2 


36 


— 






30.58 


no-fit 


3 


208 


0.66 


22.01±0.16 


2.07±0.56 


31.58 


fitted 


4 


55 


0.94 


21.90±0.31 


2.96±2.61 


30.50 


fitted 


5 


28 


— 






30.48 


no-fit 


6 


41 


— 






30.64 


no-fit 


7 


47 


— 






31.29 


Trl6-19 


8 


54 


— 






30.76 


no-fit 


9 


45 


— 






30.68 


no-fit 


10 


25 


— 






30.42 


no-fit 


11 


41 


— 






30.72 


no-fit 


12 


132 


1.50 


21.23±0.62 


2.65±0.94 


31.08 


fitted 


13 


89 


0.40 


21.64±0.14 


0.67+0.93 


30.83 


fitted 


14 


15 


— 






30.23 


no-fit 


15 


41 


— 






30.65 


no-fit 


16 


23 


— 






30.40 


no-fit 


17 


28 


— 






30.48 


no-fit 


18 


31 


0.80 


22.32±0.32 


1.36±0.88 


31.04 


fitted 


19 


7 


— 






29.93 


no-fit 


20 


30 


— 






30.51 


no-fit 


21 


25 


1.21 


21.97±0.65 


0.29±0.18 


30.02 


fitted 


22 


176 


0.94 


21.81±0.22 


3.32±1.45 


31.46 


fitted 


23 


25 


— 







30.43 


no-fit 


24 


22 


— 






30.38 


no-fit 


25 


177 


1.87 


22.00±0.18 


2.17±0.74 


31.52 


fitted 


26 


44 


0.90 


21.70±0.27 


1.49±0.54 


30.38 


fitted 


27 


51 


0.50 


22.06±0.32 


2.11±1.27 


31.03 


fitted 


28 


77 


0.47 


21.28±0.72 


1.49+0.44 


30.83 


fitted 


29 


36 


— 






30.59 


no-fit 


30 


12 


— 






30.12 


no-fit 



Notes: last column flag: Sources with no spectral information (no-fit) 
have their X-ray luminosities computed by using an average count-rate 
to Lx conversion factor (see §|7]l. "Hard tail" flag refers to the need of 
more components in the spectral models. Spectral fit parameters for 
OB-type stars are not listed here as they are presented in Table|4] 



atoms cm ^ mag ' . We also test the rela tion Ay =0.56 Nh+0.23 



^'* Seven sources (#150, #305, #342, #388, #692, #868 and #944) 
appear to have Nh below 10^' cm"^, and are likely foreground stars. 



[Nh in 10^'] dPredehl & Schmittlll995b . for which the median 
Avis 3.23 mag. Both of these values are in good agreement 
with the median Ay =3.6 mag computed from our near-IR anal- 
ysis (see § 4.3). The Icr dispersion of the Nh distribution is 0.4 
dex. It is translated into typical Ay range between 1.3 to 6.7 
mag. of visual extinction. X-ray sources without near-IR coun- 
terparts seem to be distributed towards higher absorption values 
(median log(NH)~21.9 cm"^) with respect to those with near- 
IR counterparts. Unfortunately, they generally have poor X-ray 
photon statistics and consequently a less reliable estimation of 
their X-ray spectral parameters. 

Unlike the Nh distribution, the kT distribution of plasma 
temperatures is not log-normal. It peaks at ~ 1.6 keV, has a me- 
dian ~ 1 .95 keV and shows an extended hard tail attributed both 
to variable sources (logPKs < -3) with harder spectra (me- 
dian kT=3.25 keV, as expected from coronal heating processes 
involved in flare-like activity) and to highly absorbed sources 
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Fig. 7. Left: The solid histogram shows the Nh (AJ^^'"''^) distribution of 415 spectral fitted sources with 2MASS counterparts. The 
upper scale was converted using an Nh/Av ratio of 1.6x10^'. The thick dashed histogram shows the A"*^"""^ mag. distribution for 
all near-IR counterparts as was computed in §|4] Both distributions peak at AJ5'*'""^=3.35 (lcr=0.4) mag and A""""^=3.26 mag. 
(lcr=0.3 dex). The shadowed histogram shows the Nh distribution of 148 sources with no near-IR counterpart. Right: Same as 
the left panel for plasma temperaturess (kT) with, in addition, the distribution for "flaring" sources (black-filled histogram). The 
peak of the overal distribution is at ~ 1 .6 keV. 



(i.e. Av >6 mag) showing a median kT~2.6 keV (while those 
in the range 1 < Ay <6 mag being softer and distributed with a 
median kT=1.75 keV). 

Finally, suspected single massive stars show typically soft 
spectra with median kT=0.62 keV, while for known massive 
binaries this value rises to kT=2.1 keV, no doubt due to hard 
X-ray photons being produced in the colliding wind region 
(CWR) of the massive O-hOB binaries. Details of X-ray spec- 
tral characteristics of massive stars in the region are presented 
in §[8] 

7. X-ray luminosity of stars 

Unabsorbed X-ray luminosities were computed for those 
sources with available spectral fits, for the [0.5-8.0] keV en- 
ergy range. For sources with no available and/or reliable spec- 
tral fit (1 19H-353 out of 1035) L^ were computed using a single 
count-rate to Lx conversion factor (CFO CF in Trumpler 16 is 
8.35°2^xlO''^ ergs/ph. Upper and lower values correspond to 
+ lcr uncertainties, respectively. 

In Figure[8] we show the X-ray luminosity distribution for 



low-maso stars, while known OB stars are plotted separately 
in the upper inset histogram. The Lx distribution of sources 



^^ It was computed as the median ratio between the individual un- 
absorbed X-ray luminosities (from best-fit spectral models) and the 
source count-rates. 

-^ We defined the low mass range M<7Mg, based on the presence 
of a significant convection envelope (at the adopted age) that explains 
X-ray emission from magnetic activity. 



has been plotted separately for: /- 592 X-ray sources with 
near-IR counterpart, not including 28 OB stars and variable 
sources (thin solid histogram), with a median Iog(Lx ) ~30.5 
ergs '. The peak of the distribution marks indeed the com- 
pleteness limit of our X-ray observation; //- 354 unidentified 
sources, except for variable ones, which appear systematically 
less luminous than those with 2MASS counterpart, with me- 
dian log(Lx ) ~30.2 erg s"' ; Hi- 11 variable sources, showing a 
median log(Lx)~3 1.1 ergs"', i.e. about 4 times higher than the 
observed Lx for similar stars in a quiescent (non-flaring) phase. 
The upper inset in Fig|8] shows the unabsorbed Lx for 28 mas- 
sive stars (including 13 binary systems) in the region. By fol- 
lowing results and discussion presented in §[8] binaries appears 
with typical LxOver lO^"'^ ergs"', higher than observed for 
those suspected single B-type stars. As we will discuss in §[8] 
massive binaries have typical X-ray luminosities Lx > 10^"*^ 
erg s"', higher than those of single OB stars. 

We examine how X-ray activity depends on stellar mass 
for Trumpler 16 low-mass stars, and compare the results with 
those already known for the ONC and CygOB2 SFRs. In 
order to increase statistics, we use Lx values computed by 
means of the average CF. While Lx and mass for Cyg OB2 
stars were computed following the same procedures used here 
(lAIbacete Colombo et aI.Ll2007l ) , fo r ONC stars they were com- 
puted diff'erently in the literatura^l For the sake of homogene- 



^' In the ONC, stellar masses were computed from optical spectra, 
available for many stars, while L^ was obtained for many stars from 
an X-ray spectral ana lysis based on a relatively high photon statistics 
JGetman et"aill2005t) . 
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Fig. 8. The X-ray luminosity distribution for 592 low-mass 
stars computed from a single Lx -Count rate conversion fac- 
tor is shown as a solid line. The grey filled histogram shows the 
Lx distribution for 354 X-ray sources without near-IR counter- 
part. The 77 variable sources are shown with the black-filled 
histogram. Upper-right inset (Hot stars) shows the unabsorbed 
Lx from X-ray spectral fits (see §[8]l. Error bars are Icr Poisson 
errors. 
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sources (Preibisch et al. . 2005 ). The thick gray line shows the 
linear regression fit to the Trl6 low-mass (0.7-2.5 M©) stars. 



and filled circles, 
tral fits, respectively. 



ity, we have re-computed Lx and masses of ONC stars by using 
a single count-rate to Lx conversion factor (CFonc=7. 52x10^-'^ 
erg/ph) and 2MASS photometry, respectively. 

In Figure |9] we show a plot of LxVS. star mass. We 
use circles for all 510 Trl6 sources with estimated masses. 
Filled (329) and open (181) circles indicate Lx values com- 
puted from spectral fits and using the CF, respectively. We 
perform a Unear regression for all sources in the 0.7-2.5 
Mo range: log(Lx)=3O.26(+O.ll)+l.O(+O.O9)log(M/M0) with 
a standard deviation in the residuals of 0.38 dex. The power- 
law slo2e_we_^iid_Jiereisin_agreement with that found 
by lAlbacete Colombo et al] (l2007h for the Cyg OB2 region: 
log(Lx)=3O.33(+O.16)+O.71(+O.13)log(M/M0) for masses in 
the 0.5-3.0 Mq range. This slope also agrees with that we 
find for the ONC stars in the mass range 0.5-3.0 M©, namely 
0.82+0.09. 

Changes in the X-ray activity of stars with differ- 
ent ages of the SFRs has been previously reported by 
Preibisch & FeigelsonI (l2005h . To address this issue, we 
computed Lx detection limits for low-mass stars in the 
Trl6(3Myrs), CygOB2 (2 Myrs) and ONC(lMyrs) observa- 
tions, as log(Lx)~30.5, 30.3 and 28.5 ergs"', respectively, 
above which a source is detectable anywhere in the FOV. 
A mass-dependent completeness fraction (fcomp) for our sur- 
vey of Trumpler 16 stars was then computed by adopti ng the 
X-ray Luminosity Functions given by Preibisch & Feigelson 



~55% in the 0.9-1.2 Mq range. Our survey of Tr 16 is statisti- 
cally complete for masses > 1.5 M©. In FigurefTOlwe present 
Lx cumulative distributions for Trumpler 16, CygOB2 and 
ONC stars in the mass range 1.5-2.5 Mq where all three data 
sets are complete. While Trumpler 16 and Cyg OB2 show very 
similar distributions, the ONC looks quite different with re- 
spect to Trl6. The distance between two distributions (D=0.31) 
was computed usi ng a two-sample Kolmogorov-Smirnov test 



(IPress et al.L 119921) . We are able to confirm, with a probability 
>99.9% (log(Pks)=-4.15), that young ONC(l Myr) stars, with 
masses in the range 1.5-2.5 Mq, are intrinsically more luminous 
in X-rays than their Trl6 (3 Myr) counterparts. 



8. X-rays from massive stars 

A variety of diff'erent physical mechanisms are responsible for 
the observed X-ray emission in OB-type stars. The most widely 
accepted explanation invokes multiple small-scale shocks in 
the inner layers of th eir radiation-driven stellar winds (e.g. 



Feldmeier et al.Lll997l) . In recent years it has gained importance 



a plasma heating m odel known as maenetically channeled 



wind shock (MCWS) (ISchulz et al.L 12003 



Owockietal..2005 



Gagneetal.Ll2005bl) . Moreover, WRh-OB and/or Oh-OB inter- 



(120051) : fcomp is -5% for 0.1-0.5 Mq, -40% for 0.5-0.9 Mq and 



acting wind binary s ystems produce an excess of X-ray emis- 
sion from the CWR ( Stevens et al.l Il992t IZhekov & Skinner , 
2000l:|Pittard & Stevenslboojl 



Table 4. Stellar parameters and X-ray spectral results for massive stars. 



# 



Name 



Stellar parameters 



Spectral 
Type 



Lbol 

[Lo] 



M 

[Mo/yr] 



[Lo] 



X-ray spectral parameters Notes 

xilA.o.f N^J Abund. kT\ kT^ Ui 

[1022]cin-2 [Zq] [keV] [keV] [erg/s] 

1.01/28 055+025 (LO) 0.63+0.17 1.95^-2 10'' 

0.98/28 075+033 0.18 0.49+0.28 2.27+1.26 ^isf^io^' 

1.9/240 O63+0.03 054 O71+0.01 2.66+0.06 5.23f^Vfi^ fl-6.7 keV 

095/54 1.50+O11 0.53 0.43+0.06 1.09+0.31 LOl^-^^lO" em-lines 

1.02/17 0.79+0.23 (1.0) 0.42+0.24 3.05+1.48 4.05f§^10'' variable 

0.90/19 O35+0.31 (1.0) 059+013 2.85+1.50 1.88^-^^10^' 

1.08/33 021+007 0.42 0.50+0.06 6.15^^*10^' 

2.12/74 0.66+0.04 1.2 0.26+0.02 1.73+0.60 1.89^^^10^^ (A),primaryCW 

" " " " " " (B),secondary,CW 

1.09/14 0.66+0.31 (1.9) 0.33+0.18 2.923-*il03' 

fiot fitted 

1.07/45 054+012 1.0 0.29+0.05 0.78+0.11 3.12f™1032 

1.08/23 058+014 (1.0) 0.31+0.07 1.29^-^1032 (A),primary 

" " " " " " (B), secondary 

(C),open 

1.13/18 023+009 (1.0) 2.87+0.88 1.26|-«103' variable 

3.75/342 7.39+039 4.25 0.27+0.32 7.16+3.6 Amf^lQ^^ (A),fl-6.7 keV 

(B),CW,variable 

0.78/23 060+009 (1.0) 0.46+0.07 5.53*-^{l03' 

1.02/13 048+029 (1.0) O30+0.16 1.0+O23 3.78f2j;i03i 

1.80/49 0.54+0.05 (1.0) 0.25+0.03 0.52+0.09 6.16^-^^1032 em-lines 

0.74/16 025+015 (1.0) 0.33+0.15 LIO^-^^IO^' 

1.08/17 093+079 (1.0) 0.11+0.17 0.96+0.74 -i.n^f^lQ^^ (A) primary 

" " " " " " (B) secondary 

1.3/64 1.07+O10 (1.9) 0.35+0.07 1.86+0.12 7.71^-^^1032 prob. binary 

not fitted 

0.93/32 0.78+0.09 (1.0) 0.28+0.03 >4.0 2.963-^2io32 gm-lines 

fiot fitted 

0.88/28 058+012 (1.0) 0.28+0.03 0.90+0.17 1.40|-^^1032 (A),SB2 

(B)SB2 

1.05/33 0.74+0.09 (1.0) 0.25+0.04 0.87+0.12 3.56^-«^1032 (A),(SB2),em-lines 

(B),SB2 

(C),SB1 

0.87/36 0.87+0.08 (1.0) 0.25+0.07 0.56+0.11 5.48^-|l032 em-lines 

1.01/10 054+032 (1.0) 0.28+0.18 0.54+0.32 3.06^|^103i 

0.83/29 1.15+015 (1.0) 0.75+0.10 >7.5 



7 

52 

74 

89 

136 

207 

228 

242 

242 

281 

286 

352 

407 

407 

407 

489 

649 

649 

687 

688 

689 

707 

729 

729 

730 

759 

803 

807 

808 

808 

812 

812 

812 

854 

888 

997 

997 



Trl6-19 

Trl6-124 

HD93162 

Trl6-244 

Trl6-ll 

Trl6-10 

HD93204 

HD93205A 

HD93205B 

Tr 16-21 

Trl6-14 

CPD-592600 

CPD-592603A 

CPD-592603B 

CPD-592603C 

Trl6-5 

T/Carinae -MS 

T/Carinae -B 

Tr 16-23 

Trl6-9 

HDE303308 

Trl6-3 

CPD-592628A 

CPD-592628B 

Tr 16-22 

Tr 16-74 

HD93343 

Tr 16-76 

CPD-592635A 

CPD-592635B 

CPD-592636A 

CPD-592636B 

CPD-592636C 

CPD-592641 

Trl6-115 

F015A 

F015B 



09.5V 

BIV 

WN6ha 

04If 

B1.5V 

BOV 

05V((f)) 

03.5V((f)) 

08V 

08V 

B05V 

06V((f)) 

07V((f)) 

09.5V 

B0.2V 

BIV 

02If 

05V 

07V 

09.5V 

04V((f)) 

08.5V 

09.5V 

B03V 

08.5V 

BIV 

07V(n) 

B2V 

08V 

09.5V 

07V 

08V 

09V 

05V 

08.5V 

05.5V((f)) 

08.5V 



47863 

23442 

1659586 

851138 
19055 
37154 

309029 

575440 
91201 
91201 
29512 

208929 

138038 

47863 

37153 

23442 

3715354 

309029 

138038 
47863 

467735 
72443 
47862 
37153 
72443 
23442 

138038 
15488 
91201 
47862 

138038 
91201 
58884 

309029 
72443 

309029 
91201 
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Notes: it "em-lines" reefers to stars with emission 
■well-separated components. Components (A), (B), 



lines in their optical spectra. "CW" notes those binaries with a colliding wind region. "Open" reefers to those long period binaries, i.e. with 
(C) refer to system components sorted according to mass decrease. 
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Fig. 10. Cumulative distributions of X-ray luminosities for 
Trumpler 16 low-mass stars with masses 1.5< M/Mq <2.5 
(thick solid line). Dotted and dashed lines represent the 
Lx cumulative distributions for ONC and CygOB2, respec- 
tively, in the same mass range. 
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Fig. 11. log(Lx /Lboi ) versus log(Lboi ) relation for massive stars 
in Trl6. Open boxes refer to suspected single stars and filled 
circles to known binaries. Notes: /- Because of flare-like vari- 
ability (see §5.1) observed in Trl6-#5 (B1.5V) and Trl6- 
#11 (BIV), their X-ray emission is probably dominated by 
a unknown low-mass star companion. //- The star Trl6-#22 
(08.5V-I-?) has the highest log(Lx/Lboi) ratio (~-5.56) among 
all MS OB-type stars of the region. 

Given the relatively large number of massive OB-type stars 
in the region, its appropiate to show in Table|4] most rele- 



vants stellar parameters of the massive stars (Lboi , M, and 
Lwind computed as ^Mv^O We computed Ljby use of an 
absorbed (wabs) thermal plasma model (apec). Metal abun- 
dance was fixed at Z = 0.3Zo in fitting faint sources (< 100 
ph.), while it was left as a free parameter for the remaining 
cases (see Tabled. In FigfTTI we show the Lx /Lboi relation. 
In spite of the observed scatter, one sees that the median 
Lx/Lboi for binaries (filled circles) is about 7 times larger 
(-8.3x10^^) than that of suspected single stars (-1.1x10^^0 
Four of the 13 known binaries in the analysis are well sepa- 
rated systems: the LBV //Carinae , a probable long period bi- 
nary wifli P~2026+2 days (DaminelH et al., 2007), the Wolf- 
Raye t star WR25 (HD93 162), a 208-days period binary sys- 
tem JGamen et all |2006|) and HDE 303308, resolved as a bi- 
nary system with a component separation of about ~38 AU 
projected along the fine guidance sensor (FGS) y-axio The 
stars Trl6#9 (09.5 V-i-?) and Trl6#23 (07 V-n?) show photo- 
metric variability p robably related to a secondary component 
(INelan et al.ll2004) . and thus we consider them also to be bina- 
ries. 

r 



(e.g. 'Sew ard et al. 



Besides of the several works 
1979; Seward & Chlebowski, 1982; Corcoran etal 



Albacete Colombo et al 
Sanchawala et al.L 



I 



120031: lEvans et al.L 12003 . 



1995 



I2004t 



2007h about the 



origin of the observed 
X-ray emission on massive stars of this region, this goal is 
far out of the scope of our paper and not extensively discused 
here. 



9. Summary and conclusion 

We report here results of a deep Chandra X-ray obser- 
vation pointed toward the ~3Myr old star forming re- 
gion Trumpler 16. Source detection was performed using 
the PWDdetect code, identifying 1035 X-ray sources in the 
17'xl7' AClS-1 FOV. Most of these seem to be outside the ob- 
scured V-shaped region of dust and gas. Star formation in this 
part of the masked region has probably been disrupted and/or 
diminished as the stellar winds are blocked inside the cloud, 
due to the efficiency of photo-evaporation processes caused by 
interactions of the nearby hot massive stars with the dense dust 
and gas structures. 

Data extraction was performed using the semi-automated 
IDL-based ACIS Extract package, which is well suited to the 
analysis of observations of crowded fields such as ours. 

The X-ray source list was cross-identified with optical and 
near-lR (2MASS) catalogs: 28 X-ray sources (of 44 within the 
FOV) were identified with optically characterized OB members 
of Trumpler 16 and 760 with 2MASS sources. Among these lat- 
ter sources almost all are believed to be Trumpler 16 members. 
About 90 X-ray sources without optical/NlR counterparts are 



^* Adopted values were taken from lSmitj ( l2006h -(table 1) 
^' Because of observed flare-like variability in Trl6-#5 and #11, we 
suspect that most of X-rays emission comes from a low-mass compan- 
ion, and therefore we discard these sources in computing the median 
L, . 

^^ It was unresolved along the FGS x-axis down to 20 AU, suggest- 
ing an eccentric orbit. 
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estimated to be of extragalactic nature (AGNs), while the re- 
maining X-ray sources with no counterpart are likely associ- 
ated with members that are fainter than the 2MASS complete- 
ness limit. 

In order to characterize the previously unidentified likely 
cluster members with NIR counterparts, we placed them 
on NIR color-magnitude (K, vs. H-Ks) and color-color (H- 
Ks vs. J-H) diagrams. A first estimate of interstellar extinction 
was obtained adopting a 3 Myr isochrone for the low- and 
intermediate-mass stars and assuming that O- and early B-type 
stars lie on the MS. We find a median visual absorption for OB 
stars of Av ~2.0 mag, while low mass likely members seem 
to be slightly more absorbed. Ay ~3.6 mag. We also use the 
3 Myr isochrone and the J magnitude to estimate masses of 
likely members assuming that they share the same distance and 
absorption. Our sample of X-ray selected members with near- 
IR counterparts reaches down to M=O.5-O.6M0, and is likely 
complete down to ~ 1.5 Mq. From the H-Ks vs. J-H diagram we 
estimate that ~I5% (51/339) of low-mass stars have NIR ex- 
cesses, finding it to be quite a high percentage with respect 
to the 2 Myr old Cyg OB2. We believe that the disk fraction 
in young SFRs is more dependent on spatial morphology of 
gas and dust around massive stars, which may enhance photo- 
evaporation, and thus shorten disk lifetimes, than on the total 
number of massive stars in the region. 

At least 77 sources, i.e. ~7A%, were found to be vari- 
able within our observation with a confidence level greater than 
99.9%. Only three of the 28 detected O- and early B-type stars 
were detected as variable during our 90-ksec observation, in 
spite of the high statistics of the OB stars' light curves. These 
exceptions are the known binary 03.5V-I-08V star HD93205 
(our source #242), showing a rather linear decay of the count 
rate during the observation plus a short-term variability, the 
BI.5 V star Trl6-ll (source #136) and the BIV star Trl6-5 
(source#489). The latter two show a flare-like variability prob- 
ably related to unresolved low-mass companions. 

We modeled the ACIS X-ray spectra of sources with more 
than 20 photons and where fcont <1- We assumed an ab- 
sorbed single-component thermal emission model. The me- 
dian log(NH) of the sources is 21.73 (cm~^). This value agrees 
well with the median Ay computed from the near-IR diagram. 
The median kT of low-mass stars is 2.6 keV. Sources associ- 
ated with O- and early B-type stars are instead quite soft (me- 
dian kT: 0.60 keV). Absorption corrected X-ray luminosities 
of OB stars were calculated from the best-fit spectral models. 
O and B-type stars are the most luminous, with Lx= 2.5x10-^"- 
6.3 xlO-'-' erg s"'. Their X-ray and bolometric luminosities are 
in rough agreement with the relation Lx/Lboi = 10"^, albeit 
with an order of magnitude dispersion. Low mass stars have 
Lx ranging between 10-'*' and 10-^' ergs"'(median Lx = 2.8 x 
10^°). Variable low mass stars are on average 0.5 dex brighter 
(Iog(Lx)~31.0 ergs '). These X-ray luminosities are consis- 
tent with those of similar mass stars in the slightly younger 
(2 Myr) Cyg OB2 region. However, in the mass range 1.5-2.5 
Mq, the ONC (I Myr) shows higher X-ray activity level than 
observed in Trumpler 16 stars in the same mass range. We be- 
lieve that the age-Lx activity connection is an acceptable expla- 
nation of this result. 
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